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ABSTRACT 


The  XR-3  Loads  and  Motions  Program,  a  digital  computer 
simulation  of  the  six  degrees-of- freedom  equations  of  motion 
for  the  XR-3  captive  air  iDubble  test  craft,  is  subjected  to  a 
wide  range  of  regular  and  irregular  wave  excitations.   The 
operating  characteristics  of  the  simulated  craft  and  the 
frequency  response  functions  for  the  heave,  pitch,  and  roll 
responses  in  ahead,  abeam,  and  following  seas  are  obtained 
for  regular  wave  excitations.   The  ahead  seas  frequency  response 
functions  are  verified  with  irregular  wave  excitation. 
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I.   INTRODUCTION 
A .   BACKGROUND 

Since  1973  work  has  been  ongoing  at  the  Naval  Postgraduate 
School  on  a  digital  computer  model  that  simulates  the  forces 
generated  in  and  the  motion  experienced  by  the  school's  two 
ton  XR-3  captive  air  bubble  test  craft.   The  origin  of  this 
model  was  the  "SES  Loads  and  Motions  Program"  developed  for 
the  U.S.  Navy  by  Oceanics  Inc.  for  the  Bell  100  ton  (100-B) 
surface  effect  ship  /~Ref .  1_7.   Since  there  were  major  design 
differences  between  the  100-B  and  the  XR-3,  Leo  and  Boncal 
/~Ref .  2_7  converted  the  100-B  model  to  an  initial  representa- 
tion of  the  XR-3  test  craft. 

In  197^  Finley  /~Ref .  3_7  refined  the  bow  and  stern  seal 
subroutines  of  the  XR-3  model  and  developed  new  lift  fan  maps 
in  order  to  improve  the  model's  performance  by  providing  more 
realistic  air  flow  rates.   At  the  same  time  Forbes  /~Ref.  hj 
measured  various  variables  of  the  actual  XR-3  test  craft  while 
in  operation  and  compared  them  with  the  computed  values  from 
the  XR-3  model.   Forbes  also  introduced  several  modifications 
to  various  subroutines  of  the  model.   By  the  end  of  197^  there 
were  three  versions  of  a  program  to  simulate  the  XR-3. 

During  1975  Menzel  /~Ref .  Sj   undertook  a  comparison  of 
Finley* s  and  Forbes 's  versions  by  looking  at  their  pitch  and 
roll  characteristics.   Based  on  the  results  of  this  study, 
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Menzel  introduced  program  modifications  to  improve  the  com- 
puted pitch  and  roll  behavior  of  the  simulated  XR-3.   Menzel* s 
version  is  considered  to  be  the  best  to  date  but  there  are 
still  portions  of  the  model  that  require  study. 

Before  additional  modifications  are  made,  however,  it  was 
deemed  desirable  to  provide  a  baseline  with  which  further 
versions  could  be  compared.   In  addition  it  was  noted  that 
the  majority  of  the  past  work  had  been  performed  while  oper- 
ating the  model  in  calm  water  conditions.   Although  the  model 
includes  a  "sea-state  generator,"  relatively  few  runs  had  been 
made  simulating  the  XR-3  craft  in  the  presence  of  waves.   Two 
reasons  account  for  this.   Past  work  has  been  primarily  con- 
cerned with  developing  the  force  interactions  within  the 
model  where  the  response  of  the  model  to  wave  forcing  actions 
was  not  desired.   Second,  a  large  amount  of  CPU  time  was 
necessary  to  simulate  a  few  seconds  of  simulated  run  time. 
The  response  of  the  craft,  however,  to  various  wave  inputs 
provides  a  commonly  used  baseline  for  future  comparisons. 

B.   OBJECTIVE 

The  objective  of  this  study  was  to  investigate  the  dynamic 
and  operational  characteristics  of  Menzel* s  model  when  sub- 
ject to  many  differing  operating  conditions.   In  particular 
it  was  desired  to  obtain  the  frequency  response  of  the  model, 
i.e.,  the  response  of  such  characteristics  as  heave,  pitch, 
and  roll  as  functions  of  wave  frequency  and  amplitude.   In 
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support  of  this  primary  objective  it  was  desired  to  review 
the  characteristics  of  ocean  waves  and  the  interaction  of 
waves  and  craft  in  the  XR-3  digital  model. 
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II.   FREQUENCY  RESPONSE  FUNCTIONS 

The  fundamental  principle  of  linear  systems  is  that  in 
the  frequency  domain  the  response  of  a  system,  Y(w)  ,  is  equal 
to  the  excitation,  X(w) ,  multiplied  "by  the  system's  transfer 
function,  H(w) : 

Y(w)   =  X(w)H(w) . 

Also       |y(w)|   =|X(w)|  |h(w)| 

/Y(w)    =  /X(w)  +  /H(w) 


where  |  |  denotes  the  magnitude  and  /__   the  phase  and 

Sy(w)  =  S^(w)-|h(w)|^ 

where  S  (w)  and  S  (w)  represent  either  the  power  spectral 
X        y 

density  or  energy  spectral  density  of  the  excitation  and 

response  respectively. 

This  concept,  ingrained  in  the  electrical  engineer,  is 

also  used  by  naval  architects  and  marine  engineers  to  describe 

the  motion  of  a  displacement  type  vessel  in  a  random  sea  and 

was  first  proposed  for  such  use  by  St.  Denis  and  Pierson 

/~Ref.  ^J'      "^^^  marine  engineer  refers  to  the  transfer  func- 

I     I  2 
tion  as  the  "frequency  response  function"  and  |  H(w)|    as 

the  "Response  Amplitude  Operator,"  (RAO).   The  excitation  is 

the  ocean. wave  and  the  excitation  frequency  is  given  in  terms 

of  the  encounter  frequency  (see  Appendix  A).   Using  this 


16 


concept,  the  vessel  can  be  represented  as  a  single  input, 
multiple  output  system  as  shown  in  Figure  1.   The  frequency 
response  functions  thus  describe  the  vessel  in  terms  of  sev- 
eral linear  input-output  relationships.   There  can  be  a  sep- 
arate frequency  response  function  for  each  degree  of  freedom. 
They  can  be  obtained  either  theoretically  or  experimentally 
by  such  techniques  as  towing  a  model  in  regular  waves  in  a 
towing  tank  and  measuring  the  response  of  interest. 


R-L(t) 


n(t) 


wave 
displacement 


>  Responses 


Ri(t) 


RN^t)  ^ 


Frequency 

Response 

Functions 

Craft  Input-Output  Relationship 


Figure  1 
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The  assumption  that  the  excitation-response  relationship 
is  linear  is  not  true  for  the  SES  craft  since  the  hydrodynamic 
and  aerostatic  characteristics  of  the  SES  are  inherently  non- 
linear.  However  a  theoretical  frequency  response  model  for 
the  XR-3  was  developed  by  the  Aerojet-General  Corp.  /~Ref .  ?_/ 
using  a  four  degree-of- freedom  representation  of  the  craft 
(responses  in  surge,  heave,  pitch,  and  roll).   It  was  assumed 
that  the  response  motions  were  small  enough  so  that  the  devia- 
tion from  the  linear  response  was  negligible  for  the  motions 
in  the  vertical  plane. 

The  frequency  response  characteristics  obtained  from  this 
model  are  reproduced  in  Appendix  B  and  discussed  further  in 
Chapter  IV.   The  predictions  were  compared  to  the  results  from 
i   scale  model  towing  tests  and,  in  general,  the  experimental 
and  theoretical  frequency  response  information  was  found  to 
show  good  agreement. 

The  purpose  of  this  study  was  to  obtain  the  frequency 
response  functions  for  the  XR-3  Loads  and  Motions  program,  a 
time-domain  model  that  simulates  the  craft  in  six  degrees-of- 
freedom  and  provides  time  domain  outputs  of  the  craft's  motions 
and  forces  for  the  on-bubble  mode  of  operation.   Since  the 
model  provides  time-domain  outputs,  determination  of  the 
response  characteristics  was  basically  an  experimental  process 
consisting  of  running  the  model  in  a  sinusoidal  wave  many 
times  for  different  wave  frequencies,  amplitudes,  and  craft 
speeds.   The  desired  responses  were  observed,  measured,  and 
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plotted  as  functions  of  the  encounter  frequency.   The  process 
and  results  are  discussed  in  detail  in  Chapter  IV. 

However,  because  the  program  provides  for  all  six  degrees- 
of-freedom,  the  inherent  nonlinearities  of  the  SES  craft  are 
incorporated  into  the  model,  particularly  in  the  coupling  of 
the  various  degrees  of  freedom.   The  generality  of  the  fre- 
quency response  obtained  is  therefore  dependant  on  the  signi- 
ficance of  the  nonlinearities. 
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III.   WAVE  SIMULATION  IN  THE  MODEL 

The  XR-3  Loads  and  Motions  program  has  the  capability  of 
simulating  the  craft's  motion  in  regular  and  irregular  seas. 
The  purpose  of  this  section  is  to  review  "briefly  the  charac- 
teristics of  ocean  waves  and  their  modeling  in  the  computer 
program. 

A.   WAVE  CHARACTERISTICS 

The  outstanding  characteristic  of  the  open  ocean  is  its 
irregularity  which  can  he  described  hy  statistical  methods 
on  the  assumption  that  a  large  number  of  regular  (sinusoidal) 
"waves  having  different  wavelengths,  frequencies,  directions, 
phase,  and  amplitudes  are  superimposed  to  form  the  randomly 
varying  sea.   Consequently  a  brief  look  at  the  properties  of 
regular  waves  is  necessary. 

It  is  common  to  characterize  wave  motion  by  the  direction 
of  the  field  displacement;  waves  are  either  transverse  or 
longitudinal.   In  the  transverse  wave  the  field  intensity 
oscillates  in  a  direction  perpendicular  to  the  direction  of 
wave  travel,  as  in  the  electromagnetic  wave.   In  the  longi- 
tudinal wave  the  particles  of  the  medium  move  in  a  direction 
parallel  to  the  motion  of  the  wave  itself  as  in  the  sound 
wave.   However  ocean  waves  can  be  considered  a  combination 
of  the  motion  experienced  in  both  the  transverse  and  longi- 
tudinal waves  since  the  path  of  an  individual  water  particle 
is  circular  as  a  wave  moves  by. 
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The  properties  of  regular  waves  are  well  defined  by 
gravity-free  surface  wave  theories.   It  is  assumed  that  the 
wave  crests  are  straight,  infinitely  long,  parallel  and 
equally  spaced  and  that  wave  heights  are  constant.   The  wave 
advances  in  a  direction  perpendicular  to  the  line  of  the 
crests  at  a  uniforTii  velocity,  C,  often  referred  to  as  "cel- 
erity" to  emphasize  that  it  is  the  wave  forTn  rather  than  the 
water  particles  that  advances. 

For  convenience  the  regular  wave  is  normally  modeled  as 
a  sinusoid  with  the  following  properties: 

Wave  Velocity:   C  =  (-f^)*  =  -f^  =  |-        (IH-l) 
Frequency:      w  =  (^^)*  =  ^^  (III- 2) 


Wavelength:     ^  =   ^  ^^^  (III-3) 

w 

Wave  surface  profile: 

A(x,t)  =  A^  cos  (^-J^  -  wt  +  £.  )     (III-4) 

O  A 

where   w  =  angular  frequency  (radians/sec) 

9<  =  wavelength 

T  =  wave  period 

g  =  gravitational  constant 

X  =  distance  along  the  direction  of 
propagation  from  a  defined  origin 

£.  =  phase  angle  at  the  origin 

A  =  maximum  wave  amplitude 
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The  slope  of  the  wave  surface  is  obtained  by  differentiating 
(III-'^): 


A(x,t)    ^   "^^^o   <,,.,/2jZLX 


sin(^^^-^  -  wt  +^) 

^  A^  A 

The  maximum  wave  slope  is  then: 

S-ttA  360   A 

V  =  — radians  =  — degrees  (III-5) 

A  A 

An  important  property  of  the  regular  wave  is  the  energy 
associated  with  it  consisting  of  kinetic  energy  associated 
with  the  orbital  motion  of  the  water  particles  and  potential 
energy  resulting  from  the  change  in  water  levels  in  wave 
hallows  and  crests.   The  average  energy  per  unit  area  of 
surface  for  a  regular  wave  is  /~Ref .  8_7: 

2 

E  =  i  k  A^  (III-6) 

where  k  is  a  water  density  constant. 

The  irregular  elevation  or  amplitude  of  the  sea  surface  is 
considered  to  be  composed  of  a  theoretically  infinite  number 
of  small  amplitude  sine  waves  of  different  frequencies.   Since 
the  wave  frequencies  differ,  the  velocities  of  propagation 
differ  (III-l)  so  that  at  a  particular  spot  on  the  ocean  the 
phases  of  the  individual  waves  are  uniformly  distributed  over 
2  7t  radians.   The  wave  amplitude  at  that  spot  is  the  sum  of 
the  individual  wave  amplitudes  and  is  therefore  also  random. 

Oceanographers  describe  the  random  amplitude  of  the  sea 
by  use  of  an  energy  density  spectrum.   This  indicates  the 
amount  of  energy  in  the  different  component  waves  which  have 
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combined  to  produce  the  observed  irregular  pattern.   It  is 
a  plot  of  the  energy  density  per  unit  frequency,  S(w),  versus 
frequency,  (w) .   S(w)  is  equal  to  E/  6 w  where  E  is  the  energy 
per  unit  area  of  the  wave  surface  and  <i  w  is  an  incremental 
frequency  element. 

The  area  under  the  spectrum  represents  the  total  energy 
per  unit  surface  area  of  the  sea  surface  since: 

.w 


dw  =  /   S(w)dw  =  E 
o  ^  o. 

The  area  between  any  two  frequencies,  w, ,  and  Wp,  represent 
the  energy  density  associated  with  all  waves  between  the  two 
frequencies.   It  can  be  shown  that  the  total  energy  in  this 
frequency  increment  ^  w  at  a  central  frequency  w  is  /~Ref . 

kS(w  )  <6w. 

Since  the  regular  waves  the  energy  per  unit  surface  area  is 
proportional  to  the  square  of  the  wave  amplitude  (III-6)  it 
follows  that  the  square  of  the  amplitude  of  a  wave  having  the 
same  energy  as  all  the  wave  components  in  a  band  of  frequencies 
represented  by  ,^  w  is: 

A^  =  2S(w^)  6^  (III-7) 

This  allows  the  construction  of  an  approximate  energy  density 
spectrum  from  a  finite  number  of  regular  wave  components  for 
a  given  sea  condition. 
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Certain  visible  properties  of  the  random  sea  that  can  be 
measured  "by  trained  observers  have  also  been  related  to  the 
total  energy  E  and  are  listed  in  Table  I  /~Ref .  9j,      The 
most  commonly  referred  to  property  is  the  Significant  Wave 
Height,  defined  as  the  average  of  the  heights  of  the  highest 
third  of  the  waves  observed  and  usually  written  as  H-.  z^. 

Several  mathematical  models  of  the  energy  density  spectrum 
have  been  proposed  but  the  most  widely  used  model  today  for 
wind  generated  sea  states  is  the  Pierson-Moskowitz  spectrum 
/~Ref .  10_7'   It  is  given  by  the  equation 

S(w)  =  Q-QQ^^  ^  ^ ^  (III-8) 

w-^ 

.  where  w^  =  (0  .  210 -g/H-j^ /^)^. 

and  is  plotted  in  Figure  2  as  a  function  of  significant  wave 
height. 

Pierson  and  Moskowitz  also  developed  a  relationship 
between  the  significant  wave  height  and  wind  velocity,  U, 
as: 

H^/3  =  0.210  ^.  (III-9) 

TABLE   I 
Approximate  Wave  Height  Data 

Most   frequent  wave  height    1 . ^IjE 

Average  wave  height    1 • 77/E 

H-L/3    2.83/e 

H    /        3.6o/e" 
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Figure  2,   Pierson-Moskowitz  Energy  Density  Spectrums 
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Substitution  of  equation  III-9  into  III-8  yields 

S(w)  =  2^221^ -^ (III-IO) 

It  is  this  form  of  the  Pierson-Moskowitz  spectrum  that  can 
be  used  in  the  XR-3  model  to  generate  wave  components  for 
the  simulation  of  a  random  sea  as  discussed  in  the  follow- 
ing section  and  Chapter  V. 

B.   SIMULATION  IN  THE  MODEL 

Two  subroutines  of  the  XR-3  Loads  and  Motions  program, 
INCON  and  WAVES,  are  involved  in  the  construction  and  simu- 
lation of  waves.   INCON  develops  the  wave  components  to  be 
used  in  the  construction  of  the  regular  or  irregular  seas 
and  prepares  the  wave  parameter  table  section  of  the  program's 
summary  listing.  WAVES  calculates  the  wave  forces  and  moments 
acting  on  the  craft  and  generates  the  wave  amplitudes  at  the 
various  stations  around  the  hull  as  well  as  the  bubble  volume 
lost  due  to  wave  elevation. 

1.   Subroutine  INCON 

Four  options  are  available  to  define  the  wave  compon- 
ent frequencies  and  amplitudes. 

Option  1:   The  amplitudes  and  frequencies  (rad/sec) 

of  up  to  10  wave  components  are  specified. 

Option  2:   The  amplitudes  and  wavelengths  of  up  to 
10  wave  components  are  specified. 

Option  3'      The  significant  wave  height  (H-,  /o)  and 

significant  range  of  periods  are-^specif  ied . 
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Option  ^:   The  significant  wave  height  and  range  of 

significant  wave  frequencies  are  specified. 

Options  1  and  2  allow  the  user  to  construct  the  regular  wave 
with  one  component  or  irregular  seas  composed  of  up  to  10 
components.   In  option  3  and  k   the  program  calculates  the 
wave  components  from  the  Pierson-Moskowitz  energy  density 
spectrum  based  on  the  H-,  /o  and  the  range  of  frequencies  in 
the  spectrum.   If  option  3  is  used  the  upper  and  lower  period 
limits  specified  are  converted  to  the  lower  (w,)  and  upper 
(w  )  frequency  limits  of  the  spectrum. 

The  desired  number  of  component  amplitudes  and  fre- 
quencies are  developed  from  the  spectrum  as  follows.   The 
frequency  range  is  divided  into  as  many  increments  as  there 
are  desired  wave  components.   The  width  of  each  component  is 
w, c^  where  n  is  the  increment  number  and  c  is  a  weighted 
scaling  factor  defined  as 

w  "-/^ 

where  N  is  the  number  of  wave  components. 

The  central  frequency  of  each  increment,  w^,  is  deter- 
mined by  finding  the  midpoint  of  the  increment , and  the  energy 
density  per  unit  frequency,  S(w) ,  is  calculated  from  equation 
III-IO  for  each  w  .  The  central  frequencies  of  each  of  the 
increments  are  then  used  as  the  component  frequencies.  The 
component  wave  amplitudes  are  found  by  taking  the  square  root 
of  equation  III-7.  i.e., 

A^  =  (2S(w^) -h   w)2. 
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The  above  computations  are  performed  in  INCON  from  statement 
INCN36OO  to  INCN^IOO. 

After  the  frequency  and  amplitude  of  each  wave  compo- 
nent is  determined,  the  waveslope,  initial  encounter  frequency, 
and  encounter  period  of  each  wave  component  is  calculated 
(statements  INCN^720  through  INCN^920).   This  information  is 
listed  in  the  wave  parameters  tahle  section  of  the  program's 
summary  listing.   It  should  "be  stressed  that  the  encounter 
frequencies  listed  are  only  true  if  during  the  simulation  the 
speed  is  maintained  at  the  initial  condition  speed. 
2.   Subroutine  WAVES 

Once  the  regular  wave  components  with  the  appropriate 
distribution  of  amplitude  and  frequency  have  been  determined 
in  INCON,  subroutine  WAVES  constructs  the  irregular  sea  simply 
by  adding  together  the  regular  wave  components  as  time  pro- 
gresses.  It's  important  to  note  that  this  is  not  a  true 
random  sea  since  the  irregular  wave  formed  by  the  sum  of 
several  periodic  regular  waves  must  also  be  periodic  unless 
of  course  one  of  the  regular  wave  frequencies  is  an  irrational 
number.   However,  it  is  also  apparent  that  the  period  of  the 
irregular  wave  may  be  exceptionally  long  as  compared  to  that 
of  any  of  the  regular  waves  and  that  in  the  course  of  the 
simulation,  one  cycle  of  the  irregular  wave  may  not  be  com- 
pleted.  Therefore  the  irregular  sea  developed  by  adding 
together  the  component  regular  waves  may  be  considered  pseudo- 
random. 
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WAVES  generates  the  wave  amplitudes  at  the  various 
stations  around  the  sidewalls  and  seals  relative  to  calm 
water.   The  wave  elevation  relative  to  any  point  on  the  craft 
is  represented  for  the  n   component  wave  by 

ETA^(x,y)  =  A^  sin  ^   (-xcosT  -  y  sinT  )  +  f^(t)  (Ill-ll) 


^n 


where 

f^it)    =  w^t  +  ^  (-X^  cos/9    -  Y^   sin/3    )      (III-12) 

Xq  and  Y  represent  the  position  of  the  center  of 
gravity  of  the  craft  with  respect  to  the  inertial  reference 
frame  and  are  time-vaiying  quantities  as  the  craft  moves. 
/3  is  the  angle  between  the  normal  to  the  wave  crests  and  the 
inertial  frame  X  axis  as  shown  in  Figure  3-   The  craft  axes 
are  x  and  y  and  the  two  define  any  point  on  the  craft  in  the 
horizontal  plane.   The  angle  Y  is  the  angle  between  the  ship 
X  axis  and  the  normal  to  the  wave  crests  and  is  defined  as 
y  =  /3  -  '&' 

It  is  not  difficult  to  see  how  the  relative  motion 
between  the  craft  and  wave  is  accounted  for  in  equation  III-ll. 
Assume  that  the  wave  height  at  the  center  of  gravity  is  desired 
for  a  craft  heading  straight  into  the  waves.   Then  x  =  0, 
y=  0,  /3=   180°,  and  "€    =   0° .      Equations  III-ll  and  III-12 
then  reduce  to  one  equation 

2-7r 


ETA^(c.g.)  =  A^  sin  (w^t  +  ±^   x^)  . 


29 


wave  crest 


/ 

wave  direction 


Figure  3.   Relation  between  ship  axes,  wave 
direction,  and  inertial  reference 
frame. 
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But  Xq  is  Vt,  where  Y  is  the  craft  velocity,  and  w  =  ^^ 

n 
(Eq.  III-2)  giving 

ETA^(c.g.)  =  A^  sin(  ^^   +  -^)t. 
n         "       ^n     ^n 

The  argument  of  the  sine  is  therefore  the  encounter  frequency 

times  time  were  w   =  (C  +  V)  ^^  as  shown  in  equation  A-3. 

^n  ^ 

Equations  III-ll  and  III-12  are  represented  by  statements 

WAVS2530,  WAVS05^0,  and  WAVSO9IO . 

Because  the  integration  routine  used  in  the  program 

cannot  calculate  an  instantaneous  rise  in  the  "build-up  of 

the  wave  amplitude,  the  amplitude  of  the  wave  is  restrained 

by  a  factor 

AMPFAC  =  1.0  -  e-^/'^-^^^S?^ 

This  restriction  of  the  wave  amplitude  is  evident  in  the  time 
responses  shown  in  chapter  k. 

The  primary  function  of  WAVES  is  to  calculate  the 
wave  forces  and  moments  acting  on  the  craft.   This  is  described 
in  detail  in  reference  1.   WAVES  also  has  the  option  of 
printing  the  wave  elevations  at  the  various  stations  and 
center  of  gravity,  the  plenum  volume  lost  due  to  the  waves  , 
and  the  wave  forces  and  moments. 
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IV.   REGULAR  SEA  RESPONSE  CHARACTERISTICS 
A.   PROCEDURE 

Determination  of  the  response  characteristics  of  the  XR-3 
Loads  and  Motions  computer  program  was  basically  an  experimen- 
tal process.   In  brief,  a  particular  value  of  craft  speed,  wave 
amplitude,  and  wave  heading  was  chosen  and  held  constant  while 
the  model  was  repeatedly  run  in  a  sinusoidal  regular  wave, 
each  run  having  a  new  wave  frequency.   The  desired  oscillatory 
response  was  measured  for  each  run   and  plotted  as  a  function 
of  the  encounter  frequency.   A  new  value  of  either  craft  veloc- 
ity, wave  amplitude,  or  wave  direction  was  chosen  and  the  runs 
were  repeated. 

The  craft  velocities  used  were  10,  20,  and  30  knots  which 
covered  the  speed  range  of  the  actual  XR-3  test  craft.   Three 
wave  headings  of  180,  90,  and  0  degrees  were  chosen  thereby 
simulating  the  craft  in  ahead,  abeam,  and  following  seas 
respectively.   The  wave  frequencies  were  chosen  to  cover  the 
range  of  frequencies  contained  in  a  random  sea  and  present 
integer  values  of  encounter  frequency  to  the  craft  as  dis- 
cussed in  Appendix  A.   The  responses  of  primary  interest 
were  roll,  pitch,  and  heave  at  the  center  of  gravity. 

The  XR-3  model  can  be  run  with  constant  speed  or  constant 
thrust  options.  With  the  constamt  speed  option  the  thrust  is 
varied  as  necessary  to  maintain  the  required  speed.  With  the 
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constant  thrust  option,  the  thrust  is  held  constant  and  the 
speed  is  allowed  to  vary  as  the  simulation  progresses.   The 
data  presented  here  was  obtained  while  operating  the  program 
with  the  constant  thrust  option  as  it  simulates  the  operation 
of  the  actual  XR-3  test  craft  more  closely.   However  the  con- 
stant speed  option  was  used  frequently  to  obtain  initial 
operating  conditions  as  is  explained  later  in  this  chapter. 

The  program's  initial  conditions  in  pitch,  draft,  plenum 
pressure,  and  thrust  were  determined  as  a  function  of  speed 
and  are  given  in  Appendix  C.   Appendix  D  lists  the  craft 
parameters  used.   The  model  was  normally  run  for  35  seconds 
to  ensure  that  the  steady-state  responses  had  been  obtained. 
The  responses  were  assumed  to  have  reached  steady-state  when 
their  peak-to-peak  and  average  values  remained  constant  for 
at  least  ten  seconds.   The  responses  were  normally  sampled 
and  printed  every  0.05  seconds  to  ensure  that  the  minimum 
Nyquist  rate  was  more  than  satisfied  at  the  higher  encounter 
frequencies. 

1.   Ahead  Seas  Runs 

The  initial  runs  were  performed  with  a  craft  speed  of 
20  knots  in  ahead  seas  {  (B  =   180°).   A  wave  amplitude  of  0.2 
feet  was  chosen  as  the  excitation  and  the  20  knot  calm  water 
thrust  of  186.72  lbs.  per  side  was  used.   Steady-state  values 
of  heave  and  pitch  were  obtained  rapidly  for  encounter  fre- 
quencies of  1,  2,  and  3  radians  per  second.   However  at 
w  =  ^  steady-state  responses  could  not  be  obtained  merely 


33 


by  running  the  model  for  a  greater  period  of  time.   As  shown 
in  Figures  ^-A  through  ^-C,  as  time  progressed  the  draft 
increased  significantly  while  the  pitch  oscillations  diverged. 
Concurrently  the  speed  had  dropped  in  60  seconds  from  an  ini- 
tial 20  knots  to  approximately  15.6  knots  thereby  reducing 
the  frequency  of  encounter  with  the  wave  excitation.   More 
important  was  the  fact  that  the  plenum  volume  and  pressure 
were  found  to  have  decreased  significantly  as  shown  in  Figures 
^-D  and  ^-E. 

Intuitively  the  above  can  be  explained  as  follows.   As 
the  plenum  volume  and  pressure  fell,  the  craft  sank  deeper 
into  the  water,  increasing  the  drag  on  the  sidewalls  and  reduc- 
ing the  speed.   The  only  possible  source  for  the  loss  of  pres- 
sure and  volume  was  believed  to  be  associated  with  the  leakage 
of  air  through  gaps  formed  between  the  sidewalls  and  the 
water  surface  or  the  bow  seal  and  water  surface.   It  was 
found  that  the  most  forward  section  of  the  sidewall,  and 
therefore  presumably  the  bow  seal  also,  did  rise  in  and  out 
of  the  water  allowing  air  leakage  through  the  gap  formed  as 
discussed  in  Appendix  E. 

However  it  should  be  noted  that  the  pressure-volume 
relationship  inside  the  plenum  is  an  adiabatic  process  which 
would  have  the  pressure  increase  instead  of  decrease  for 
constant  air  mass  with  a  reduced  volume.   Therefore  for 
pressure  and  volume  to  both  drop,  the  mass  of  the  air  must 
be  decreasing  at  a  rate  faster  than  the  volume.   Since  the 
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Figure  k-B,      Pitch  vs.  Time 
V=20  kts,   A=0.2  ft,   w^=^ 
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Figure  ^-D.   Plenum  Pressure  vs.  Time 
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Figure  ^-E.   Plenum  Volume  vs.  Time 
V=20  kts,   A=0.2  ft,   Wg=^ 
Constant  Thrust,   /3  =180° 
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air  flow  rate  into  the  plenum  will  saturate  due  to  the  air 
fan  map  limitations,  high  leakage  rates  would  cause  the  mass 
of  the  air  to  drop  faster  than  the  volume  and  thereby  cause 
the  pressure  to  drop. 

Of  interest  also  is  the  initial  transient  portion  of 
the  responses.   As  the  wave  amplitude  rises  exponentially, 
the  responses  are  almost  linear.   However  as  the  coupled  heave 
and  pitch  amplitudes  increase,  gapping  "begins,  resulting  in 
the  rapid  drop  in  volume  and  therefore  causing  the  nonlinear 
saturation  of  the  air  fans  which  is  reflected  in  the  pressure 
response. 

Clearly  the  calm  water  thmst  value  was  not  sufficient 
to  overcome  the  increased  drag  due  to  the  increasing  draft. 
Because  of  this,  the  model  was  ran  with  the  constant  speed 
option  in  order  to  determine  the  thrust  necessary  to  main- 
tain the  20  knot  speed  and  the  encounter  frequency  of  ^  radians 
per  second.   The  result  is  shown  in  Figure  5  where  the  nec- 
essary thrust  is  a  not  quite  sinusoidal  variable  with  an 
average  value  of  215  lbs.  and  lagging  the  wave  excitation  by 
approximately  2^0  degrees.   The  model  was  then  run  with  a 
constant  thrust  value  of  215  lbs.  per  side  and  steady-state 
responses  were  obtained  as  shown  in  Figure  6.   Note  that  the 
craft  still  sinks  somewhat  deeper  into  the  water  than  for 
calm  water  conditions  and  that  the  plenum  pressure  and  volume 
have  reached  a  steady-state  condition. 
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Figure  6-A.   Heave  vs.  Time 
V=20  kts,   A=0.2  ft,  ^^=^' 
Constant  Thrust=215  lb,  /3   =180° 
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Figure  6-B.   Pitch  vs.  Time 
V=20  kts,   A=0.2  ft,   Wg=^ 
Constant  Thrust  =215  1^,  ^   =l8o' 
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Figure  6-C.   Speed  vs.  Time 
V=20  kts,   A=0.2  ft,   Wg=^ 
Constant  Thrust  =215  11:^,  /3  =l8o' 
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From  the  initial  20  knot  runs  a  general  procedure  was 
developed  for  obtaining  steady-state  responses.   The  model  was 
initially  run  with  the  constant  thrust  option  using  the  calm 
water  thrust  for  the  desired  speed.   The  pitch,  heave,  and 
roll  responses  and  craft  speed  were  observed.   If  steady- 
state  values  were  not  reached  within  thirty-five  seconds  the 
model  was  rerun  with  the  constant  speed  option.   The  average 
value  of  thrust  from  this  run  was  then  used  as  the  thrust  for 
a  second  constant  thrust  run.   This  would  normally  produce 
steady-state  heave  and  pitch  responses  for  speeds  of  20  and 
30  knots. 

Obtaining  steady-state  responses  for  a  craft  speed  of 
10  knots  required  several  more  iterations  of  the  thrust,  how- 
ever.  For  example,  at  an  encounter  frequency  of  six  radians 
per  second,  the  calm  water  thrust  of  200.31  lbs.  per  side 
drove  the  craft  to  a  speed  of  11.5  knots  in  20  seconds  and 
resulted  in  a  decreasing  draft  and  converging  heave  and  pitch 
oscillations.   A  constant  speed  run  produced  an  average  thrust 
of  19^  lbs.  per  side  but  when  this  was  applied  as  the  input 
to  a  constant  thrust  run,  the  speed  dropped  below  hump  speed. 
A  choice  midway  between  the  two  values  (I98  lbs.)  again  pro- 
duced an  increasing  speed  and  unsteady  oscillations.   A  final 
thrust  value  of  I96  lbs.  produced  the  desired  steady-state 
responses. 

All  runs  performed  at  10  knots  were  found  to  be 
extremely  sensitive  to  small  changes  in  thrust.   This  is  due 
to  operating  the  program  in  a  region  where  the  drag  versus 
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speed  has  a  minima,  allowing  small  variations  in  thrust  to 
significantly  effect  the  craft  speed.   This  problem  was  also 
noted  by  Cagle  /f"Ref .  llJ7  in  his  study  of  the  original  100-B 
digital  model. 

2.   Abeam  Seas  Runs 

To  obtain  the  frequency  response  functions  for  the 
roll  characteristics  of  the  digital  model,  the  program  was 
run  simulating  the  craft  in  abeam  seas,  i.e.,  the  angle  bet- 
ween the  normal  to  the  wave  crests  (in  the  direction  of  pro- 
pagation) and  the  craft's  X  -  axis  was  90  degrees.   The  craft 
was  run  at  20  and  30  knots  using  the  calm  water  thrust  values 
and  the  constant  thrust  option.   The  wave  amplitude  was  0.2 
feet. 

A  typical  set  of  responses  obtained  while  operating 
the  program  in  this  manner  is  shown  in  Figures  7-A  through 
7-E  where  neither  the  heave,  pitch,  nor  roll  reached  steady- 
state.   Figure  7-E  provided  the  explanation  for  the  responses. 
As  the  craft  encountered  the  waves  a  yaw  oscillation  was 
induced.   As  time  progressed  the  average  value  of  the  yaw 
became  non-zero  and  the  craft  gradually  turned  into  the  wave, 
changing  the  wave  excitation-craft  response  relationship. 
To  obtain  steady-state  pitch,  roll,  and  heave  conditions, 
it  was  clear  that  it  was  necessary  to  first  have  a  stable 
yaw  response. 

To  obtain  steady-state  yaw  a  rudder  map,  i.e.,  a  set 
of  data  points  that  describe  the  rudder  angle  as  a  function 
of  time,  was  added  to  the  model  in  order  to  turn  the  craft 
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Figure  7-C.  Roll  vs.  Time 
V=20  kts,  A=0.2  ft,  Wg=3 
/3  =  90°,   Constant  Thrust 
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Figure  7-D.   Speed  vs.  Time 
V=20  kts,   A=0.2  ft,   Wg=3 
x9  =  90°,   Constant  Thrust 
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in  a  direction  opposite  that  which  was  induced  by  the  wave. 
Several  iterations  of  this  process  were  required  before  steady- 
state  responses  were  reached.   It  was  determined  that  a  small  • 
rudder  angle,  typically  less  than  a  degree,  of  constant  mag- 
nitude would  allow  steady-state  operation  to  be  reached.   For 
the  example  given  in  Figure  7,  a  constant  rudder  map  of  -0.2 
degrees  as  shown  in  Figure  8-A  produced  steady-state  responses 
as  exemplified  in  Figure  8-B  and  8-C.   Although  the  average 
value  of  the  yaw  is  still  not  zero,  its  magnitude  is  consider- 
ably less  than  initially  when  the  rudder  angle  was  zero. 
3.   Following  Seas  Runs 

Following  seas  are  encountered  when  the  craft  and  waves 
travel  in  the  same  direction.   The  procedures  described  for 
ahead  seas  were  used  to  determine  the  desired  response  char- 
acteristics. 

B .   RESULTS 

The  responses  of  primary  interest  were  heave,  pitch,  and 
roll.   Once  steady-state  conditions  had  been  reached,  the 
peak-to-peak  and  average  values  were  determined  and  are  listed 
in  Appendix  G.   Since  heave  is  a  function  of  vertical  wave 
displacement,  the  ratio  of  the  two  was  obtained  and  plotted 
as  a  function  of  encounter  frequency.   Similiarily  the  ratio 
of  pitch  amplitude  and  roll  amplitude  to  their  excitation, 
the  wave  slope,  was  found  and  plotted. 


5^ 


ft 

a 

m 

u 

<D 

TJ 

TJ 

13 

cn 

c8 

C^ 

0) 

II 

■^ 

S 

<U 

CO 

•H 

^ 

^ 

EH 

• 

" 

EH 

to 

-P 

> 

ch 

-P 

^ 

CM 

cfl 

cd 

• 

+J 

>H 

o 

w 

il 

i:: 

< 

o 

• 

o 

< 

1 

1 

00 

w 

■p 

o 

o 

^ 

o 

?-l 

ON 

a 

o 

C\J 

II 

•H 

II 

/VN 

P4 

> 

^ 

(se9j5ap)      /wex 


55 


<D     C^ 


B 

II 

-P 

•H 

<U 

CQ 

EH 

^ 

3 

• 

^ 

CQ 

•• 

EH 

> 

4^ 

ch 

•P 

H 

^ 

H 

OJ 

a 

O 

• 

p 

« 

o 

CO 

11 

s:^ 

< 

o 

• 

o 

pq 

CO 

m 

•P 

o 

<D 

^ 

o 

u 

ON 

a 

o 

CM 

II 

•H 

II 

m 

fo 

> 

^ 

(sesjtSep)      lio^ 


56 


0) 

S 

•  H 

Eh 

CO 

> 

r-\ 

O 

;h 

-p 

c 

o 

o 

^ 

0) 

n 

cn 

T3 

II 

•P 

^ 

OJ 

m 

« 

^ 

^ 
^ 

^ 

x: 

-p 

« 

en 

•H 

-p 

^ 

l+H 

-p 

5 

CM 

aJ 

CtJ 

• 

-P 

>H 

o 

m 

II 

^ 

< 

o 

• 

o 

o 

1 

^ 

OD 

w 

•> 

-p 

o 

<U 

^ 

o 

U 

ON 

13 

o 

tU) 

CM 

11 

•H 

II 
> 

^ 

(seejcSap)      avbx 


57 


1.   Frequency  Response  Functions 
a.   Ahead  Seas 

Program  rims  simulating  ahead  seas  conditions 
over  a  range  of  encounter  frequencies  from  1  to  20  radians 
per  second  were  conducted  at  the  following  craft  speeds  and 
wave  heights. 

Speed  Wave  Amplitude 

10  knots  0.2  ft 

20  knots  0.1,  0.2,  0.3,  0.5  ft  ' 

30  knots  0.2,  0.3  ft 

Given  that  with  a  sinusoidal  excitation  a  linear  system  should 
produce  a  sinusoidal  response,  in  general  the  heave  and  pitch 
responses  were  found  to  be  linear  for  the  ahead  seas  condi- 
tions.  Small  nonlinearities  were  noted  at  w  =1  in  the  pitch 
response,  two  examples  of  which  are  shown  in  Figure  9-      The 
only  extreme  occurrence  was  at  V=20  knots,  w  =4,  and  A=0.5  ft. 
As  shown  in  Figure  10  "both  pitch  and  heave  displayed  large 
nonlinear  response  characteristics. 

Figures  11  through  1?  present  the  frequency 
response  functions  for  heave  and  pitch  for  the  seven  combi- 
nations of  craft  speed  and  wave  heights  simulated.   The 
predominate  characteristic  in  all  of  them  is  the  resonant 
peak  in  the  pitch  response  near  w  =^.   A  comparison  of  the 
responses  as  a  function  of  waveheight  for  speeds  of  20  and 
30  knots  is  shown  in  Figures  18  and  19.   The  heave  response 
is  nearly  independent  of  waveheight  except  for  the  region 
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Figure  11. 
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V=10  kts,    A=0.2  ft, 

Ahead  Seas 

Heave  Frequency  Response 

Pitch  Frequency  Response 
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Figure  12.   V=20  kts,   A=0.1  ft, 
Ahead  Seas 
Q  =  Heave  Frequency  Response 
0  =  Pitch  Frequency  Response 
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Figure  I3.   V=20  kts,   A=0.2  ft, 
Ahead  Seas, 
CD  =  Heave  Frequency  Response 
0  =  Pitch  Frequency  Response 
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Figure  1^.   V=20  kts,   A=0.3  ft, 
Ahead  Seas 
Q  =  Heave  Frequency  Response 
O  =  Pitch  Frequency  Response 
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Figure  l6.   V=30  kts,   A=0.2  ft, 
Ahead  Seas 
Q  =  Heave  Frequency  Response 
O  =  Pitch  Frequency  Response 
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Figure  17.   V=30  kts,   A=0.3  ft, 
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Q  =  Heave  Frequency  Response 
0  =  Pitch  Frequency  Response 
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Figure  I8-B.   Pitch  Frequency  Response  vs.  Wave 
Height,  Y=20  kts 
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Figure  I9-A.   Heave  Frequency  Response  vs.  Wave 
Height,  Y=30  kts 
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Figure  19 -B.   Pitch  Frequency  Response  vs.  Wave 
Height,  V=30  kts 
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near  Wg=^.   This  is  not  true  for  the  pitch  response  where  at 
the  higher  frequencies  the  response  falls  off  faster  for 
greater  wave  heights.  A   comparison  of  the  responses  as  a 
function  of  speed  for  constant  waveheight  is  shown  in  Fig- 
ure 20.   Clearly  the  response  in  both  heave  and  pitch  is 
dependent  on  the  speed  of  the  craft. 

A  comparison  with  the  four  degree-of- freedom 
theoretical  response  in  Figure  B-1  in  Appendix  B  is  also 
of  interest.   Both  models  show  the  speed  dependence  of  heave 
but  not  of  pitch.   In  general  the  theoretical  response  char- 
acteristics fall  off  slower  in  both  pitch  and  heave  then  that 
obtained  from  the  XR-3  digital  simulation  and  the  theoretical 
responses  would  appear  to  be  independent  of  waveheight  at 
resonance.   Also  noteworthy  is  the  second  peak  in  heave  near 
w  =8  for  the  four  degree-of- freedom  model  which  does  not  appear 
in  the  six  degree-of- freedom  model  until  a  much  higher  fre- 
quency is  encountered. 

The  frequency  response  functions  obtained  from 
the  XR-3  Loads  and  Motions  program  show  that  at  low  frequen- 
cies representing  waves  with  large  wavelengths  as  compared 
to  the  craft  length,  or  the  length  of  the  bubble  chamber  for 
the  CAB  type  craft,  the  simulated  XR-3  will  follow  the  wave 
surface  and  slope.   At  high  frequencies  the  wavelength  is 
small  relative  to  the  plenum  length  and  the  craft  will  tend 
to  not  follow  the  waves  but  will  "platform"  over  them.   Inter- 
mediate to  these  two  conditions  is  a  range  of  resonant  fre- 
quencies for  which  wave  inputs,  particularly  the  wave  slope, 
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Figure  20 -A.   Heave  Frequency  Response  Vs.  Speed, 
A=0.2  f-t 
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Figure  20-B 


Pitch  Frequency  Response  vs.  Speed, 
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would  "be  amplified  to  result  in  craft  motions  that  could  "be 
significantly  larger  than  the  input  values.   Figure  21  plots 
the  dependence  of  heave  and  pitch  on  the  wavelength  to  plenum 
length  ratio.   Note  that  in  terms  of  this  ratio,  heave  is 
fairly  independent  of  speed, 
"b.   Following  Seas 

Following  seas  runs  were  simulated  for  craft  speeds 
of  10  and  20  knots  with  a  wave  amplitude  of  0.2  feet.   Steady- 
state  conditions  could  not  be  obtained  for  V=10  knots,  w  =8. 

e 

Nonlinearities  were  noted  primarily  at  10  knots  in  heave  for 
the  lower  encounter  frequencies.   Figure  22  presents  repre- 
sentative examples. 

The  frequency  response  functions  for  the  following 
seas  simulations  are  presented  in  Figures  23  and  2^.   The  10 
knot  response  characteristics  can  be  compared  with  the  four  ■ 
degree-of- freedom  theoretical  response  in  Figure  B-2.   As  can 
be  seen  the  two  bear  little  resemblance  to  each  other.   In 
particular  the  magnitudes  are  considerably  smaller  for  the 
digital  simulation.   The  digital  simulation's  pitch  response 
has  the  general  features  of  the  theoretical  response,  in 
particular  the  notch  above  w  =2  is  evident  in  both,  but  the 
heave  response  of  the  digital  simulation  has  a  resonant  peak 
in  contrast  to  the  notch  in  the  theoretical  response, 
c.   Abeam  Seas 

Runs  simulating  abeam  seas  conditions  were  conducted 
with  a  wave  amplitude  of  0.2  feet  for  craft  speeds  of  20  and 
30  knots.   For  abeam  seas  the  frequency  of  encounter  is  not 
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V  =  10  kts,   A=0.2ft,   w  =1 

e 


V  =  10  kts,   A  =  0.2  ft,.  Wg  =  2 


V  =  10  kts,   A=0.2ft,   w  =^ 

e 


Figure  22.   Heave  Nonlinearities 
Following  Seas 
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Figure  23.   V=10  kts,    A=0 . 2  ft, 

Following  Seas 
Q  =  Heave  Frequency  Response 

O  =  Pitch  Frequency  Response 
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Figure  2k.      V=20  kts ,    A=0 . 2  ft 
Following  Seas 
Q  =  Heave  Frequency  Response 
O  =  Pitch  Frequency  Response 


20 


80 


a  function  of  the  craft  speed  but  is  the  actual  wave  frequency. 
This  has  the  effect  of  producing  considerably  greater  wave 
slopes  than  were  encountered  in  the  ahead  or  following  seas 
conditions. 

Nonlinearities  were  found  in  pitch  and  roll  at 
the  lowest  encounter  frequencies.   However  since  pitch  is 
coupled  to  both  the  roll  and  yaw,  nonlinearities  in  its 
response  were  to  be  expected.   Typical  responses  are  presented 
in  Figure  25- 

The  heave  frequency  responses  for  the  two  speeds 
were  so  nearly  identical  that  they  can  be  plotted  as  one  as 
in  Figure  26.   The  frequency  response  functions  for  pitch 
and  roll  are  given  in  Figures  27  and  28  and  also  show  a  high 
degree  of  independence  from  speed.   Comparisons  with  the 
theoretical  responses  are  not  possible  since  they  are  based 
on  a  craft  speed  of  10  knots.   However  the  theoretical  pitch 
response  does  deserve  some  attention.   The  pitch  frequency 
response  function  shows  that  for  the  lower  encounter  frequen- 
cies the  pitch  amplitude  is  equal  to  the  wave  slope.   If 
there  were  no  yaw,  as  believed  to  be  the  case  for  the  four 
degree-of- freedom  model,  the  pitch  of  the  craft  would  be 
minimal,  much  as  is  characterized  by  the  responses  from  the 
digital  simulation. 

2.   Operating  Characteristics 

The  operating  characteristics  of  the  XR-3  program 
were  defined  as  the  average  draft,  average  pitch,  and  the 
thrust  necessary  to  maintain  the  desired  speed.   Figures  29 
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Figure  2?.   V=20  kts,   A=0.2  ft. 
Abeam  Seas 
O  =  Pitch  Frequency  Response 
^  =  Roll  Frequency  Response 
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Figure  28.   V=30  kts,   A=0.2  ft, 
Abeam  Seas 
O  =  Pitch  Frequency  Response 
A  =  Roll  Frequency  Response 
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through  36  present  this  information  as  a  function  of  encounter 
frequency,  craft  speed,  wave  height,  and  craft  heading  with 
respect  to  the  waves.   For  both  ahead  and  following  seas  the 
20  and  30  knot  runs  show  great  similarity.   At  resonance  there 
is  a  large  increase  in  average  pitch  and  draft  along  with  a 
significant  increase  in  the  thrust  necessary  to  counter  the 
increased  drag  and  maintain  the  desired  speed.   However  the 
10  knot  runs  are  markedly  different  in  that  during  resonance 
it  was  necessary  to  decrease  the  thrust  to  maintain  the  speed. 

For  the  abeam  seas  conditions,  the  calm  water  thrust 
was  adequate  to  maintain  the  desired  craft  speed.   In  all 
cases  except  one,  the  average  draft  and  pitch  remained  at 
their  calm  water  values.   As  noted  in  the  procedures  section 
steady-state  yaw  was  a  prerequisite  for  obtaining  steady- 
state  responses  in  heave,  pitch,  and  roll.   The  peak-to-peak 
and  average  values  of  yaw  at  the  end  of  30  seconds  run  time 
without  rudder  control  as  well  as  the  peak-to-peak  value 
with  rudder  control  are  listed  in  Tables  G-10  and  G-11  of 
Appendix  G. 
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Figure  29.   Operating  Characteristics 
V=10  kts,   A=0.2  ft, 
Ahead  Seas 
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V.   IRREGULAR  SEA  RESPONSE  CHARACTERISTICS 

A  vessel's  behavior  in  regular  waves,  as  was  obtained  for 
the  digital  XR-3  in  chapter  IV,  is  of  fundamental  importance, 
for  if  the  vessel  can  be  described  by  linear  responses,  its 
motion  in  an  irregular  seaway  can  be  described  as  the  super- 
position of  the  responses  of  the  vessel  to  all  the  wave  com- 
ponents of  the  seaway,  i.e., 

Y(w)  =  X(w)  •  H(w) 

To  check  the  validity  of  the  ahead  sea  frequency  response 
characteristics  obtained  for  the  XR-3  digital  model  in  regu- 
lar seas,  the  program's  ability  to  simulate  the  craft's  motion 
in  irregular  seas  was  used. 

The  irregular  seas  were  constructed  using  wave  option  3 
in  which  given  the  H-,  /^  and  the  range  of  significant  periods 
the  component  wave  frequencies  and  amplitudes  are  calculated 
by  the  program  using  the  Pierson-Moskowitz  spectrum.  Each 
irregular  sea  generated  was  composed  of  ten  wave  components. 

The  significant  wave  heights  and  range  of  significant 
wave  periods  for  various  sea  conditions  were  obtained  from 
data  compiled  by  the  David  Taylor  Model  Basin,  a  portion  of 
which  is  listed  in  Table  II.   The  associated  spectrums  were 
plotted  in  Figure  2.  Using  the  procedures  developed  for  ahead 
sea  runs  in  regular  waves,  simulations  were  conducted  for  the 
conditions  listed  in  Table  III.   It  was  found  that  in  sea 
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TABLE  II 
Significant  Wave  Heights  and  Wave  Periods 


Sea 
State 

Wing 
Velocity 
(knots) 

^1/3 
(ft.) 

^1/10 
(ft.) 

Significant 
Range  of 
Periods  (sec) 

0 

2 

0.08 

0.10 

0-1.2 

1 

5 

0.29 

0.37 

0.'^-2.8 

1 

8.5 

1.0 

1.2 

0.8-5.0 

2 

10 

1.4 

1.8 

1.0-6.0 

2 

12 

2.2 

2.8 

1.0-7.0 

2 

13.5 

2.9 

3.7 

1.^^-7.6 

TABLE  III 
Irregular  Sea  Simulation  Conditions 


Craft 

Speed 

(knots) 


^1/3 
(ft.) 


Thrust 
per  side 
(lb) 


20 


30 


0.29 

1.0 

1.4 

1.0 

1.4 


186.72 

210.0 

223.0 

300.0 

310.0 
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states  with  significant  wave  heights  greater  than  those  listed 
in  Table  III  there  was  nearly  continuous  contact  with  the  top 
of  the  plenum  chamber,  a  condition  which  is  not  accounted  for 
in  the  program  and  therefore  avoided.   A  typical  example  of 
the^  generated  irregular  wave  as  seen  at  the  craft's  center  of 
gravity  and  the  responses  to  the  wave  are  shown  in  Figure  3?-A 
through  37-E.   Clearly  apparent  in  the  heave  and  pitch  responses 
is  the  "low  pass  filtering"  of  the  wave  excitation's  high 
frequency  components  which  was  the  dominate  characteristic 
of  the  heave  and  pitch  regular  sea  frequency  response  functions. 
The  magnitude  of  the  frequency  response  function  is 

H(w)  =  -&• 

In  order  to  determine  the  excitation  and  response  spectrums, 
the  discrete  fourier  transform  (DFT)  was  taken  of  the  wave 
excitation  using  an  IBM  36O  library  subroutine  as  described 
in  Appendix  F. 

To  obtain  an  accurate  representation  of  the  spectrums 
with  the  DFT,  several  conditions  had  to  be  met.   To  prevent 
aliasing,  the  high-frequency  components  of  the  time  function 
impersonating  low  frequency  components,  the  sampling  rate  of 
the  time  response  or  excitation  had  to  be  greater  than  the 
Nyquist  rate.   A  sampling  interval  of  O.03  seconds  satisfied 
this  condition.   The  frequency  resolution  of  the  DFT  is  l/r 
where  n:    is  the  sample  length.   The  longer  the  sample  the 
greater  the  resolution.   However  this  conflicted  with  the 
actual  running  of  the  program  in  which  greater  than  2.5  minutes 
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Figure  37-A.   C.G.  Wave  Height  vs.  Time 
V  =  20  kts,    H^ /^  =  1.0  ft 
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Figure  37-B.   Heave  vs.  Time 
V  =  20  kts,   H^  /  =  1.0  ft 
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Figure  37-E.   Plenum  Pressure  vs.  Time 
V  =  20  kts,  K^/^  =   1.0  ft 


10  2 


CHJ  time  was  required  to  simulate  the  craft  for  one  second  of 
run  time  in  irregular  seas.   Finally,  the  bandwidth  of  the 
spectrum  is  a  function  of  the  number  of  sample  points  and 
the  sample  length  /"Ref.  12_7.   A  512  point  DFT,  which  there- 
fore produced  256  frequency  components,  and  a  sample  length 
of  approximately  16  seconds  were  found  to  produce  good  discrete 
spectrums . 

Figure  38-A  shows  the  discrete  magnitude  spectrum  of  the 
wave  height  at  the  center  of  gravity  for  V=20  kts,  H-,  /^=1.0ft. 
Note  that  the  spectrum  produced  is  a  function  of  encounter 
frequency  since  the  data  was  taken  while  the  craft  was  in 
motion.   Also  shown  in  the  figure  is  the  encounter  frequencies 
and  amplitudes  of  the  components  used  in  actually  constructing 
the  irregular  wave.   The  two  agree  quite  well. 

Figures  38-B  and  38-C  present  the  magnitude  spectrum  of 

the  heave  and  pitch  responses  to  the  above  wave  excitation, 

again  as  a  function  of  encounter  frequency.   The  "low  pass" 

characteristics  of  the  heave  and  pitch  regular  sea  frequency 

response  functions  are  clearly  evident  in  the  spectrums  also. 

However  also  noted  are  the  components  at  frequencies  below 

w  =3  for  which  there  is  no  excitation.   These  values  were 
e 

ignored.   From  the  wave  spectra  data,  the  frequencies  and 
amplitudes  of  the  peaks  were  obtained.   The  amplitudes  of  the 
heave  spectra  were  also  obtained  at  the  same  frequencies  and 
the  ratio  of  the  two  were  calculated  (see  Table  G-12,  G-13. 
and  G-1^)  and  are  plotted  in  Figure  38-D.   As  a  comparison 
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Figure  38-D.   Heave  Frequency  Response 
V=20  kts,   flw^=1.0  ft 
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Figure  38-E.   Pitch  Frequency  Response 
V=20  kts,   H-^/^=1.0  ft 
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with  the  heave  frequency  response  resulting  from  regular  seas, 
the  response  for  V=20,  A=0.3  ft.  is  also  plotted. 

To  determine  the  irregular  sea  pitch  frequency  response 
it  was  necessary  to  first  calculate  the  actual  wave  frequency 
from  the  encounter  frequency  and  then  the  maximum  wave  slope 
for  each  wave  frequency.   The  ratio  of  the  pitch  amplitude  to 
the  maximum  wave  slope  was  then  determined  and  is  plotted  in 
Figure  38-E.   Note  that  for  "both  the  heave  and  pitch  irregular 
sea  frequency  response  functions  there  is  close  agreement 
with  those  obtained  in  the  regular  seas  runs. 

The  spectrums  of  wave  amplitude  and  pitch  and  heave 
responses  for  two  other  combinations  of  sea  state  and  craft 
speed  are  shown  in  Figures  39  through  ^0.   Also  presented 
are  the  associated  frequency  response  functions.   As  can  be 
seen,  the  frequency  response  functions  for  the  irregular  sea 
runs  agree  well  with  the  regular  sea  responses. 


109 


10,0 


20,0 


30.0 


^0.0 


w   (rad/sec) 


Figure  39-A.   C.G.  Wave  Height  Spectrum 
V  =  20  kts,  K^y^  =  lA   ft 
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Figure  39-B.   Heave  Spectrum 


V  =  20  kts,   H-^ /^  =  1.^  ft 
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Figure  39-D.   Heave  Frequency  Response 
V=20  kts,   H^ /^=1.^  ft 
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Figure  39-E.   Pitch  Frequency  Response 
V=20  kts,  R^y^=lA   ft 
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Figure   ^0-A.      C.G.    Wave  Height   Spectrum 
V  =   30  kts,      R^y^  =   1.0   ft 
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V  =  30  kts,   H^  /  =  1.0  ft 


116 


0.7    - 


0.6    - 


0.5   - 


CO 
<D 

a> 

u 
(no 

-p 


A 


O 
-P 
•H 


0.3    - 


0.2     - 


Ool 


10.0 


w  (r ad/sec) 


Figure  ^0-C.   Pitch  Spectrum 
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Figure  ^0-D.   Heave  Frequency  Response 
V=30  kts,   H-L/^=1.0  ft 
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Figure  ^0-E.   Pitch  Frequency  Response 
V=30  kts,   H-^  /  =1.0  ft 
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VI .   CONCLUSIONS  AND  RECOMMENDATT ONS 


The  XR-3  Loads  and  Motions  Program  was  subjected  to  a 
wide  variety  of  wave  conditions.   The  conditions  necessary 
to  obtain  steady-state  operation  were  determined  and  the 
frequency  response  functions  for  heave,  pitch,  and  roll  were 
developed.   The  heave  and  pitch  frequency  response  functions 
obtained  from  the  regular  ahead  seas  runs  were  found  to  com- 
pare well  with  those  obtained  from  irregular  ahead  seas  by 
use  of  the  DFT.   In  general  the  nonlinearities  inherent  in 
a  six  degree-of- freedom  model  were  found  to  be  minimal  and 
exist  at  the  lower  encounter  frequencies.   It  is  believed 
that  the  results  from  this  study  provide  a  good  baseline 
with  which  to  compare  future  modifications  to  the  model. 

A  most  interesting  study  would  be  to  obtain  frequency 
response  data  from  the  actual  XR-3  test  craft  and  compare 
it  with  that  obtained  from  the  digital  model.   It  is  recom- 
mended that  a  closer  look  be  given  to  the  model's  simulation 
of  the  craft  at  the  lower  speeds.   It  was  noted  that  during 
resonance  at  10  knots,  as  the  craft  sank  deeper  actually 
less  thrust  was  required  to  maintain  the  desired  speed  than 
at  the  other  encounter  frequencies.   There  appears  to  be  no 
plausible  explanation  for  this. 

Although  the  conditions  under  which  gapping  of  the  side- 
walls  occurred  were  determined  no  effort  was  made  to  determine 
the  actual  air  flow  rates  out  of  the  gaps  or  if  indeed  the 
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air  loss  was  great  and  rapid  enough  to  cause  the  loss  in  plenum 
pressure  noted  and  therefore  the  settling  of  the  craft  deeper 
into  the  water.   Finally,  although  not  analyzed  or  presented 
in  this  study,  data  was  obtained  and  is  available  on  the 
vertical  accelerations  of  the  center  of  gravity  for  the  major- 
ity of  the  steady-state  runs.   Analysis  of  this  data  could  be 
of  use  in  defining  the  habitability  zones  for  the  simulated 
XR-3. 
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APPENDIX  A 
Encounter  Frequency 

The  frequency  of  excitation  is  commonly  given  in  terms 
of  the  encounter  frequency,  the  frequency  actually  "seen"  hy 
the  vessel.  Due  to  the  forward  speed  of  the  craft,  the  fre- 
quency of  encounter  with  the  waves  is  increased  or  decreased 
depending  on  the  relative  direction  between  the  vessel's 
course  and  the  direction  of  propagation  of  the  surface  wave. 
The  encounter  frequency  is  given  by: 

w  =  -^T?-  (C-V  cos/3)  (rad/sec)  (A-1) 

e    A 

where   >  =  wave  length   (ft) 

C  =  wave  velocity  (ft/sec) 

V  =  craft  velocity  (ft/sec) 

/3  =  the  angle  between  the  craft's  course  and 

the  direction  of  propagation  of  the  wave. 

Given  that  the  wave  velocity,  C  =  (|-^)^,  and  the  wavelength, 

r^~  2l21^   where  w  is  the  wave  frequency,  by  substitution  into 
w 
A-1  it  can  be  shown  that; 

2   V  cos  3      ,  ( h-9\ 

w  =  w-w  • \^  ^) 

e  g 

For  a  vessel  heading  into  the  waves,  P  =  180   and 

w  =  2.2^  (C+V).  (A-3) 

s     A 


122 


For  this  study  it  was  desired  to  simulate  regular  waves 
that  produced  integer  values  of  w  "between  1  rad/sec  and  20 
rad/sec  for  a  given  value  of  vessel  velocity.   The  application 
of  the  quadratic  formula  to  equation  A-2  yields 


1  -  /1+  (^Vw.)/g 
w  =  ^- ^ ^ (A-^) 

g 


for  /3  =  180°. 
This  relationship  is  plotted  in  Figure  A-1 
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Figure  A-1.   Wave  Frequency  vs.  Encounter 
Frequency 
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APPENDIX  B 

Theoretical  Frequency  Response  Functions 

A  theoretical  frequency  response  model  for  the  XR-3  was 
developed  by  the  Aerojet-General  Corp.  /~Ref .  7_y'  using  a  four 
degree-of- freedom  representation  of  the  craft  with  responses 
in  heave,  pitch,  roll,  and  surge.   It  was  assumed  that  the 
response  motions  were  small  enough  so  that  deviation  from 
the  linear  response  was  negligible  for  the  motions  in  the 
vertical  plane.   The  requency  response  characteristics  obtained 
are  reproduced  in  Figures  B-1  thru  B-3  and  show  not  only  the 
theoretical  response  but  the  results  from  f  scale  model  tow- 
ing tests  as  well. 
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Figure  B-IA.   Heave  Frequency  Response, 
Ahead  Seas 
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Figure  B-IB.   Pitch  Frequency  Response, 
Ahead  Seas 
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Figure  B-2A.   Heave  Frequency  Response, 
Following  Seas 


V=io 


Figure  B-2B.   Pitch  Frequency  Response, 
Following  Seas 
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Figure  B-3A.   Heave  Frequency  Response, 
Oblique  Seas 


Figure  B-3B.   Pitch  Frequency  Response, 
Oblique  Seas 
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Figure  B-3C.   Roll  Frequency  Response, 
Oblique  Seas 
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APPENDIX  C 
Initial  Conditions 

Prior  to  irunning  the  computer  model  with  simulated  wave 
conditions  the  "calm  water"  initial  conditions  had  to  be 
determined.   Previous  work  by  Forbes  /~Ref .  kj/   and  Menzel 
/~Ref .  Sjy   had  shown  that  the  initial  conditions  used  in  the 
program  had  to  be  close  to  the  steady  state  operating  values. 
Failure  to  ensure  this  would  result  in  either  an  imbalance  in 
the  forces  generated  within  the  program  and  cause  the  program 
to  stop  due  to  an  improper  integration  step  size  or  transients 
that  obscured  the  response  of  the  variable  under  study. 

The  steady  state  conditions  found  by  Menzel  for  his 
Program  3  /Table  XVI,  Ref,  SJ7   were  used  as  the  initial  input 
to  the  computer  model.   The  program  was  run  using  the  constant 
speed  option  for  forty  seconds  simulation  time  at  nine  different 
speeds.   The  steady  state  value  of  pitch  angle,  draft,  plenum 
pressure,  and  thrust  were  determined  for  each  speed.   These 
values  were  then  used  as  the  initial  conditions  for  a  second 
set  of  constant  speed  runs  to  verify  that  the  steady  state 
values  produced  no  transients  when  used  as  initial  conditions. 
A  third  set  of  mns  with  the  constant  thrust  option  in  use 
was  performed  to  verify  that  the  speed  was  maintained  at  its 
desired  value. 
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The  "calm  water"  steady  state  values  determined  by  the 
above  procedure  and  used  throughout  this  study  as  initial 
conditions  are  listed  in  Table  C-1  and  plotted  in  Figure  C-1. 


TABLE  C-1 
Calm  Water  Initial  Conditions 


Speed 
(knots) 

Pitch 
Angle 
(deg) 

Draft 
(in) 

Plenum 
Pressure 
(psf) 

Thrust 
per  side 
(lb) 

10.0 

1.62 

8.17 

23.93 

200.31 

12.5 

1.11 

7.03 

2^.8^ 

174.71 

15.0 

0.8^ 

6.7^ 

2^'.Qk 

167.85 

17.5 

0.63 

e.i^i 

2^.8^- 

173.23 

20.0 

0.^8 

6.12 

'       2^1-. 84 

186.72 

22.5 

0.36 

5.87 

2^.84 

205.77 

25.0 

0.29 

5.66 

2^.84 

229.31 

27.5 

0.25 

5.^8 

24.84 

256.65 

30.0 

0.26 

5.3^ 

24.84 

287.22 
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APPENDIX  D 
Craft  Parameters 
The  following  craft  parameters  were  used  throughout  this 
study.   The  parameters  are  listed  by  Input  Blocks.   Also 
listed  are  the  associated  variable  names  as  used  in  the 
program. 

Integrator  Tolerances:   Block  1 

Number  of  Integrators  Used:   10 

Runga-Kutta  Variable  Step  Sized  Used 

-L  -L  -^  -7 

Initial  Step  Sizes:   1x10   ,  2x10   ,  1x10  ^,    1x10  ' , 

1x10"^,  1x10"^,  1x10"^.  1x10"^,  IxlO"'^,  1x10"^ 
Masses  and  Inertias:   Block  2 
Weight:   6050  lb   WEIGHT 

Longitudinal  Center  of  Gravity:   10.05  ft  forward  of 
transom,  XS 

Vertical  Center  of  Gravity:   2.5^  ft  above  keel,  ZS 

2 
I   ,  Mass  Moment  of  Inertia  About  X  axis:   2870.O  slug- ft 

XX 

AIXX 

2 

I   .  Mass  Moment  of  Inertia  About  Y  axis:   9320.0  slug- ft  , 

yy 

AIYY 

2 

I   .  Mass  Moment  of  Inertia  About  Z  axis:   10,580.0  slug- ft  , 
zz 

AIZZ 

I   .  Mass  Moment  of  Inertia  About  XZ  axis:   -2800.0  slug-ft  , 
xz 

AIXZ 
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Craft  Geometry:   Block  3 

Total  Craft  Length:  2k. 7   ft,   XLTOT 

Number  of  Stations  along  Port  and  Starboard  Sidewalls:  11 

Number  of  Stations  along  Bow  and  Stern  Seals:   5 

Sidewal:   Block  ^ 

Y  distance  from  centerline  to  sidewall:   5.37  ft,  YSW 
Average  Wetted  Length  of  Sidewall:   21.9  ft,  XLSW 
Leakage  Orifice  Coefficient  of  Sidewall:   0.7,  CFSW 
Cross-flow  Drag  Coefficient  of  Sidewall:   1.28,  CDSW 
Average  Beam  of  Sidewall:   0.5  ft,  AVBMSW 

Stern  Seal:   Block  5 

X  corrdinate  of  Stern  Seal  Hinge:   3-97  ft  forward  of 

transom,  XSSI 
Z  coordinate  of  Stern  Seal  Hinge:   1.875  ft  above  keel, 

ZSSI 
Base  leakage  area:   3-87  ft^,  ALEAK 
Seal  Leakage  Orifice  Coefficient:   0.9,  CFSS 
Length  of  Rear  Support  Cable:   1.875  ft 
Pressure  Differential  Between  Stern  Seal  and  Plenum 

Chamber:   1  psf,  DPSS 
Length  of  Leading  Edge  of  Seal:   ^.06  ft,  XLF 

Bow  Seal:   Block  6 

X  coordinate  of  Bow  Seal  Hinge:   23.^^  ft  forward  of 

transom,  XBSI 
Seal  Leakage  Orifice  Coefficient:   0.9,  CFBS 
Pressure  Differential  between  Bow  Seal  and  Plenum 
Chamber:   1  psf,  DPBS 


13^ 


2  coordinate  of  Bow  Seal  Hinge:  1.875   ft  above  keel,  ZBSI 
Length  of  Rear  Support  Cable:   I.875  ft 
Length  of  Leading  Edge  of  Seal:  k.06   ft,  XBF 
Base  Leakage  Area:   0.08  ft^,   BLEAK 
Length  of  Middle  Support  Cable:   I.875  ft 
Plenum:   Block  7 

Plenum  Length  at  Water  Surface:   20.0  ft,  XLBW 

Plenum  Width  at  Water  Surface:   10.0  ft,  XBBW 

X  coordinate  of  Pressure  Wave  Pivot  Point:   17.2  ft  forward 

of  transom,  XPMV 
Plenum  Length  at  Deck:   20.0  ft,   XL 
X  coordinate  of  Center  of  Pressure:   10.^  ft  forward  of 

transom,   XCPO 
Plenum  Average  Height:   1.915  ft,   BUBHGT 
Froude  Number:   0.556,   FNCRIT 
Propulsion  Geometry:   Block  8 

X  coordinate  of  Propeller  Center:   -1.275  ft  forward  of 

transom,   XPO 

Y  distance  from  centerline  to  propeller  center:  5 '55   ft, 

YPO 
Z  coordinate  of  propeller:   -0.60^  ft  above  keel,   ZPO 

Rudder  Geometry:   Block  9 

X  coordinate  of  Centroid  of  Rudder,  -1.125  ft  forward 
of  transom,   XRO 

Y  distance  from  Centerline  to  Rudder  Centroid:  5-55   ft, 

YR 
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2  coordinate  of  Centroid  of  Rudder:   -0.208  ft  above  keel, 
ZRO 

Rudder  Span:   1.21  ft,   RSPAN 

Rudder  Aspect  Ratio:   2.15,   RASPR 

Rudder  Area:   0.68  ft^,   RAREA 

Average  Thickness  Ratio  of  Rudder  Section:   0.16?,   RTC 
Aerodynamics :   Block  10 

Reference  Length:   20.0  ft,   XLAERO 

Reference  Width:   10.0  ft,   BEAM 
Thrust:   Block  l6 

Starboard  and  Port  Side  Thrust  Coefficients:   0.01,  STHS 
and  STHP 
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APPENDIX  E 

Gapping 

Subroutine  SIDEWL  calculates  the  forces  and  moments  acting 
on  the  craft  due  to  the  sidewalls  as  well  as  the  leakage  flow 
rates  associated  with  any  gaps  which  open  under  the  sidewalls 
due  either  to  craft  motion  or  waves.   A  print  switch  option 
allows  printing  the  sidewall  gaps,  imersion  depths,  and  total 
forces  and  moments  at  each  integration  step.   For  purposes 
of  this  study  in  which  35  second  runs  were  normally  conducted 
the  printed  output  of  this  print  option  was  excessive  since 
all  the  above  quantities  are  printed  each  integration  step. 
Since  it  was  desired  to  know  only  if  a  gap  occurred,  its 
location  along  the  sidewall,  and  the  magnitude,  a  simple 
program  addition  was  made  to  SIDEWL  as  follows. 


C 
C 


GAP  OR  WETTED  DRAFT  CALCULATION 


"»-  GAPT=0.0 
DO  10  J=l,2 


SDWL0^60 
SDWL0^70 

SDWLO^SO 


7  GAP(J,K)=-DDIN*(1.-(DSW(J,K))/PBHEAD 
■*■   GAPT=GAPT+GAP(J,K) 


GO  TO  10 
8  GAP(J,K)=0.0 
10  CONTINUE 

IF(GAPT.EQ.O.O)  GO  TO  12 
WRITE  (6,11)VAL(1),((GAP(I,J),J=1,11). 
,       11=1.1)  ,GAPT,((GAP(I,J),J=1,11.I=2,2) 
Ul  FORMAT('(T)' ,F10.5.11(F8.5),/, 
2'(G)',FL0.5.11(F8.5)) 


SDWLO63O 


SDWL06i]'0 
SDWLO650 
SDWLO66O 
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n 

n  TTPATTAn-C       AT^rn,  SDWLO67O 

C      LEAKAGE  AREA  SDWLO68O 

-*12  ALSW=0.0 

Gapping  was  found  to  occur  under  the  sidewall  station  1 
for  the  following  conditions: 

V  =  10  knots,   A  =  0.2  ft 

Vf  =    k 
e 

V  =  20  knots,   A  =  0.2  ft 

A  =  0.3  ft 

Wq='^,  5,    6,    8,  20 

A  =  0.5  ft 

Wq=  ^,  5,    6,  8,  10 

V  =  30  knots,   A  =  0.2  ft 

Wg=  ^,  5,    6,  8,  10 
A  =  0.3  ft 

Wg=  ^,  5,    6,  8,  10,  15 
Gapping  was  found  under  station  2  for  the  following  conditions: 

V  =  20  knots,   A  =  0.5  ft 

Wg=  5,    6,  8,  10 

V  =  30  knots,   A  =  0.3  ft 

Wg=  5.  6,  10,  15 
Gapping  was  noted  at  additional  stations  for  V=20  knots, 
A  =  0.5  ft.   No  gapping  was  noted  during  steady-state  operation 
in  abeam  seas,  however  it  was  noted  under  the  most  forward 
station  in  following  sea  runs  at  V=10,  w  ='!+  and  at  station 


2  at  V=20,  w  =^'. 
e 
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APPENDIX  F 

A  DFT  Program 

The  IBM  360  library  subroutine  HARM  which  perforins 
discrete  complex  fourier  transforms  on  a  complex  three 
dimensional  array,  was  used  to  take  the  DFT  of  the  wave 
excitation  and  the  heave  and  pitch  time  responses  from 
ahead  irregular  seas  in  order  to  determine  the  frequency 
response  functions.   Once  the  complex  fourier  transform 
of  the  excitation  and  responses  was  determined,  the 
magnitude  of  the  spectrum  was  computed  and  plotted  using 
the  library  subroutine  DRAWP.   The  program  is  listed  on 
the  following  pages . 
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CINENSICN  W(600),Z(6CO),F{600) 

CI^ENSICN  FV(514) 

CINENSICN  S<1012)iINV(1012),M(2) 

Clf'ENSICiN    TIME(600)fWH(600)  ,  ZD  (  600  )  ,  PA(  6C0  J 

CC^FLEX*a    A(1024,l,l) 

IMEGER'*4    ITB(12)/l2*0/ 

RE/L*4  PTB(28)/26*0.0/ 

C 
C 

C  RE/SC    IN    THE    POWER    CF    TV<C    (DETERMINES    THE    NUNBER 

C  CF    CATA    POINTS)       SAMPLE    LENGTH,    ANC    SAMPLING 

C  INTERVAL 

C 
C 

PE/C(5,5)KK,T,DT 
£    FCFNAT(I5,Flfl.4,F10.4) 

PIs2. 14159 

M=2**KK 

N2=l 

N3«l 

NH-M/2 
C 
C 

C  RE^C    IN    TIME,    WAVE    HEIGHT,    HEAVE,    AND    PITCH    CATA 

C 

c 

CC    10    K  =  1,M 

Re/C(5,6)TIME(K) ,WH(K» ,ZD(KJ,FA(K) 
6    FCPMAT(4F16.4) 
10    CCNTINUE 

M(1)=KK 

M(2)=0 

N(3)=0 
C 
C 

C  CALCULATE    THE    FUNDAMENTAL    FREQUENCY    AND    THE 

C  FPECUENCY    COMPONENTS 

C 
C 

FLNC=(2.0*PI)/T 

FV{l)-0.0 

CC  20  J=2,NH 

FV(J)=FV( J-ll+FUND 
2C  CCNTINUE 
C 
C 

C      SET  THE  CCNPLEX  A  MATRIX  TO  ZERC 
C 
C 

CC  2  11=1, Nl 

CC  2  iasl,N2 

CC  2  12^1, N3 

2  A(I1, 12,13)^(0.0, 0.0) 
C 

C 

C      TRANSFER  THE  WAVEHEIGHT  DATA  TO  THE  COMPLEX  A  MATRIX 

C 

C 

CC  3  11*1, Nl 

CC  3  12=1, N2 

CC  3  13-1, N3 

^ (I1,1,1)=WH(I1) 

3  CCNTINUE 
C 

C 

C      CCMFUTE  THE  DISCRETE  FCLRIER  TRANSFORM  CF  THE 

C      V«/VE  HEIGHT  DATA 

C 

C 


CALL  HARM(A,M, INV,S,-1,IF£RR) 


m-o 


c 
c 
c 
c 


c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


CALCULATE  THE  MAGNITUCE  OF  THE  FGURIER  CCWFCNENTS 
^NC  CCNVERT  THE  RESULTS  TO  INCHES 

CC  20  11=1, NH 
iO  W(I1)=C4BS(A(H,1,1)  )*24.0 

CLE^R  THE  A  MATRIX 

CC  4  11*1, Nl 
CC  4  12*1, N2 
CC  4  13=1, N2 
4  4 (I1,I2,I3)*(Q.0,0.0) 

CCNFUTE  THE  DFT  OF  THE  HEAVE  DATA  AND  THE  MAGNITUCES 
CF  THE  FOURIER  CGNFCNENTS 


CC  7  11  =  1, M 

CC  7  12  =  1, N2 

CC  7  13=1, N3 

^ (I1,1,1)=Z0(I1) 

CCNTINUE 

C/SLL  H/RM(A,M,INV,S,-1,IFERR) 

CC  50  I^lfNH 

Z(I)=2.0«CABS(A(I,1,1)} 


50 


THE  CC  VALUE  OF  HEAVE   (THE  AVER/SGE  CRAFT)  IS 
STCPEC  ANC  LATER  FRINTEC.   IT  IS  NCT  PLCTTEC 


Z/VG^Zd) 
Z(1)=0.C 


CCNFUTE  THE  CFT  OF  THE  PITCH  ANGLE  DATA  ANC  THE 
J-AGMTUCES  OF  THE  FOURIER  COMPONEATS 


CC  e  11=1, M 

CC  6  12=1, N2 

CC  8  12=1, N3 

A(I1,I2,I3)=(0.0,Q.O) 

CC  «  11  =  1, M 

CC  S  12=1, N2 

CC  c  12*1, N3 

A(I1,1,1)=PA(I1) 

CCNTINUE 

CALL  H/»Ff(A,M,INV,S,-l,IFERP) 

CC  55  J=1,NH 

F  (u)=2.C*CABS(A(J,l,l) ) 


STCPE  THE  CC  VALUE  CF  PITCH  FOR 
IS  NOT  PLOTTED 


LATER  OUTPUT.   IT 


P/VG=P(1) 
P(1)=0.Q 

WFITE(e,60) 


T,DT,M,FUND,ZAVG,PAVG 


60  FCFNAT(/,«  SAMPLE  TIME  =  ',F7.4,'  SECONDS'//,'  SANF 
6LE  INTERVAL  =  ',F6.4,'  SECGNCS'//,'  NUMBER  CF  SANFL 

FUNCANENTAL  FREQUENCY  =  ',F6.4,'  RPS' 


', 14//,' 


fiES    =       . 

£//,'    ZAVG    =    •,F7.4//,'     FAVG    =    »,F7.4) 


FLCT  THE  SPECTRUMS 
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ITE(12)=1 

nE{3)  =  5 

RTE(1)=2C«C 

P7e(2)=0.S 

C/LL  CRAWP  (NH,FV,W,ITetRTB) 

PTE(1)=10.0 

P7e(2)=C,2 

CALL  CRAWP  (NH,FV,Z,ITe,RTB) 

RTE(1)=:5.0 

RT£(2)=C.l 

C4LL  CRAWP  (NH,FV,Ff ITBiRTBI 

STCF 


1^2 


APPENDIX  G 

Data  Tables 

The  frequency  response  and  operating  characteristics  data 

presented  graphically  in  the  text  are  listed  in  Tables  G-1 

through  G-1^.   The  undefined  headings  and  symbols  are: 

Ratio  1    the  ratio  of  the  peak-to-peak  heave  to  the 
peak-to-peak  wave  height 

Ratio  2    the  ratio  of  the  peak-to-peak  pitch  amplitude 
to  twice  the  maximum  wave  slope 

Ratio  3    the  ratio  of  the  peak-to-peak  roll  amplitude 
to  twice  the  maximum  wave  slope 

A        wavelength  (ft) 

L        plenum  chamber  length,  20  ft. 

*        water  contact  with  the  top  of  the  plenum 

chamber  occurred  during  steady-state  operation 

2  calm  water  draft  (in) 

0  calm  water  pitch  (degrees) 

T  calm  water  thrust  used  per  side  (lb) 

V  craft  velocity  (knots) 

A  wave  amplitude  (ft) 

/3  angle  between  the  craft  and  wave  headings 
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